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The key notion is that of null phase curves in quantum mechanical ray and 
Hilbert spaces. Examples of such curves are developed. Our generalisation is 
shown to be essential to properly understand geometric phase results in the 
cases of coherent states and of Gaussian states. Differential geometric aspects 
of null phase curves are also briefly explored. 
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I. INTRODUCTION 



The geometric phase was originally discovered in the context of cyclic adiabatic quantum 
mechanical evolution, governed by the time- dependent Schrodinger equation with a hermi- 
tian Hamiltonian operator |IL0]. Subsequent work has shown that many of these restrictions 
can be removed. Thus the geometric phase can be defined in nonadiabatic PJ, noncyclic 
and even nonunitary evolution H. Generalization to the nonabelian case has also been 
achieved |J. Finally the kinematic approach M demonstrated that even the Schrodinger 
equation and a Hamiltonian operator are not needed for defining the geometric phase. The 
intimate relationship between geometric phase and Hamilton's theory of turns @ has also 
been brought out |§. 

An important consequence of the kinematic approach has been to bring out clearly the 
close connection between geometric phases, and a family of quantum mechanical invariants 
introduced by Bargmann || while giving a new proof of the Wigner [T(J unitary- ant iunitary 
theorem. This connection depends in an essential way upon the concept of free geodesies 
in quantum mechanical ray and Hilbert spaces, and the vanishing of geometric phases for 
these geodesies. 

The purpose of this paper is to generalise this important link between Bargmann invari- 
ants and geometric phases to the widest possible extent, by going beyond the use of free 
geodesies. The key is to characterize in a complete way those ray space curves with the 
property that the geometric phase vanishes for any connected stretch of any one of them. 
We show that this property can be translated into an elementary and elegant statement 
concerning the inner product of any two Hilbert space vectors along any lift of such a ray 
space curve. We refer to these as "null phase curves" ; and the generalisation of the familiar 
statement linking Bargmann invariants and geometric phases is achieved by replacing free 
geodesies by such curves. A free geodesic is always a null phase curve; however the latter is 
a much more general object. 

The material of this paper is arranged as follows. Section 2 recalls the basic features of the 
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kinematic approach to the geometric phase; sets up free geodesies in ray and Hilbert spaces; 
shows that the geometric phase for any free geodesic vanishes; introduces the Bargmann in- 
variants; and describes their connection to geometric phases for ray space polygons bounded 
by free geodesies. In Section 3 it is argued that it should be possible to generalise this con- 
nection. This motivates the definition and complete characterization of null phase curves 
at the Hilbert space level, the previous free geodesies being a very particular case. It is 
then shown that such curves allow us to generalise the previously stated connection to the 
widest possible extent. Section 4 defines the concept of constrained geodesies in ray and 
Hilbert spaces, the motivation being that in some situations such curves may in fact be 
null phase curves. The idea is extremely simple, namely we limit ourselves to some chosen 
submanifolds in ray (and Hilbert) space, and determine curves of minimum length lying 
within these submanifolds. Section 5 examines several interesting examples to illustrate 
these ideas: a submanifold arising out of a linear subspace of Hilbert space; coherent states 
for one degree of freedom; centered Gaussian pure states for one degree of freedom; and an 
interesting submanifold in the space of two-mode coherent states. It turns out that in the 
first case constrained geodesies are just free geodesies, while in the remaining cases they 
are very different. This shows that the generalised connection between Bargmann invari- 
ants and geometric phases enunciated in this paper is just what is needed to be physically 
interesting and appropriate. In Section 6 we present a brief discussion of these ideas in the 
differential geometric framework natural to geometric phases, and also develop a direct ray 
space description of null phase curves; while Section 7 contains concluding remarks. 



II. THE CONNECTION BASED ON FREE GEODESICS 

Let H be the Hilbert space of states of some quantum system, 1Z the associated ray 
space, and n : 7i — > 1Z the corresponding projection. We shall be dealing with (sufficiently) 
smooth parametrised curves C of unit vectors in Ti, and their images C in 1Z. A curve C is 
described as follows: 
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C = {ip(s)eH\ || ^(s) 11=1, Sl <s<s 2 }cft. (2.1) 
Its image C is a curve of pure state density matrices: 

?r[C] = C c TZ , 

C = {p(s) = ^0)V>0) f I si < s < s 2 } . (2.2) 

Any C in projecting on to a given C in TZ is a lift of the latter. In particular we have a 
horizontal lift if the vectors ip( h \s) along it are such that 

U h \s),^ h \a)\=0. (2.3) 

For any curve C C TZ a geometric phase <p g [C] is defined. Its calculation is facilitated 
by going to any lift C, calculating the total and dynamical phases for C, and taking the 
difference: 

ti[C] = C : 

(f g [C) = Cptot[C) ~ <Pdyn[C) , 

(ftot[C) = arg(^(si), ip(s 2 )) , 

<^ dyn [C] = lmjds U(s), — J . (2.4) 

In particular, if C is horizontal v^dynp] vanishes, and ^p g [C] is just y9 tot [C]. 

Now we define free geodesies in TZ and Ti. Given C in TZ and any lift C in Ti, the length 
of the former can be defined as the following nondegenerate functional: 

m -J*{ n ^r-\U.).mrr. ( ,5> 



si 



It is easy to check that the integrand here is independent of the choice of lift C; it leads 
to the well known Fubini-Study metric on TZ pTT]JT^] . Free geodesies in 7Z are those C's for 
which L[C] is a minimum for given end points. And by definition a free geodesic in Ti is any 
lift of a free geodesic in TZ. It can be shown |J that any free geodesic in TZ can be lifted to 
Ti, and the parametrisation chosen, so that it can be described as follows: 



i[)(s) = 0i cos s + 02 sin s , 

(01,0l) = (02,0 2 ) = l, (0i,0 2 )=O. (2.6) 

Thus we have here a plane two-dimensional curve determined by a pair of orthonormal vec- 
tors in 7i, an arc of a circle. It may be helpful to make the following comment concerning free 
geodesies. Given any two 'non-orthogonal' points pi, p 2 G 7Z, that is such that Tr(p!p 2 ) 7^ 0, 
we can always choose unit vectors ipi,ip 2 £ 'H projecting onto p\, p2 respectively, such that 
the inner product (ipi,ip 2 ) is real positive. Then the free geodesic (2.6) will connect ip\ and 
■02 if we take 0i = ip\ and 2 = ("02 — 0i(0i,02))/ {1 — (V'ljV^) 2 } ■ It is now clear that 
■0(0) = 0i, and 0(s) = -0 2 for s = cos _1 (-0i, 2 ) G (0,vr/2). It is clear that the curve in Ti 
given by eqn. (2.6) is horizontal; and for any two points on it with \s\ — S2 1 < vr/2, the inner 
product (0(si), ip{s 2 )) is real positive, so -0(si) and tj){s 2 ) are in phase in the Pancharatnam 
sense [TCI . From these properties of free geodesies the result B 



ip g [free geodesic in 1Z] = (2.7) 

follows. This can be exploited to connect geometric phases to Bargmann invariants. 

Let ipi, ip2i ■ ■ ■ 1 V'n be any n unit vectors in Ti., no two consecutive ones being orthogonal; 
and let pi, p 2 , ■ ■ ■ , p n be their images in 7Z. Then the corresponding n- vertex Bargmann 
invariant is defined as 

A n (0i, lf) 2 , • • • ,0>„) = C01 ,1p2) ("02, -03) • • • (lpn,lfa.) 

= Tr(p lP2 ...Pn). (2.8) 

Now we draw n free geodesies in 7Z connecting pi to p 2 , p 2 to p$, . . . , p n to p\. Thus we obtain 
an n-sided polygon in 1Z bounded by free geodesies, and we can compute the corresponding 
geometric phase. Repeatedly exploiting eqn. (2.7) we obtain the basic result ||: 
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n- vertex polygon in 7Z connecting p x to p 2 , 
Pi to pz, . . . , p n to pi by free geodesies 



= -argA n (0i,0 2 , • • • ,0n) , 
Pj = ^L j = l,2,...,n. (2.9) 
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We mention in passing that this result is of considerable conceptual as well as practical value 
0- 

In connection with the above result, the following remarks may be made. As is clear 
from equation (2.8), the phases of the individual vectors ipi, ip2 ■ ■ ■ ip n can be freely altered. 
We need only to assume that successive pairs of unit vectors are not mutually orthogonal; 
then the Bargmann invariant is non zero and has a well defined phase. 



III. A GENERALISED CONNECTION 

The definition (2.8) of the Bargmann invariant requires only the choice of the n vertices 
pi, p 2 , . . . , p n G 1Z; consecutive ones need not be connected in any way to form a closed figure. 
This suggests that the connection (2.9) between these invariants and geometric phases may 
hold more generally, not only in the case where we connect pi to p 2 , P2 to p 3 , . . . , p n to pi 
by free geodesies. We now show that this is indeed so. 

We need to characterize the most general (smooth) curves C C 7Z having the property 

<p s [any connected portion of C] = . (3.1) 

We know that if C is a free geodesic, this property does follow; but there may be (indeed there 
are) many other possibilities. We can develop a simple necessary and sufficient condition on 
C such that (3.1) holds. 

Given the curve CcK, let C^ h ' be a horizontal lift and C a general lift of C in 7i. We 
have 

C = [p{s) | p(s) t = p{s) > , p(s) 2 =p(s), Trp( S ) = l, Sl <s<s 2 }; 
C w = ^ h) (s) e H | vr (>(%)) = p(s) , ^ h \s), ^ h) (s)) = j ; 

C = e H | ip{s) = e ia{s) ip (h) (s)} . (3.2) 

Here a(s) is some (smoothly varying) phase angle. For any two points on C with parameter 
values s and s' > s we have: 



<p g \p(s) to p(s') along C] = Vtot [^ h \s) to ^ h \s') along C™] 



= arg (^ h \s),^ h \s')) 

= arg (e- iQ(s V(s), e- ia(s 'V(s')) 

= arg(^(s), ip(s')) + a(s) - a(s') . 



(3.3) 



From this result we see that the necessary and sufficient condition on C to secure the property 
(3.1) can be expressed in several equivalent ways, using either an arbitrary lift C of C or a 
horizontal lift C {h) : 

y? s [any connected portion of C] = 



Here separability is to be understood in the additive, and not in the multiplicative, sense. 
It is important to recognise that these characterizations are reparametrization invariant. Any 
curve C C 1Z obeying (3.4) will be called a "null phase curve in 7Z" ; and any lift C of such 
a C will be called a "null phase curve in H n . Free geodesies are null phase curves, but not 
necessarily conversely. 

It may be helpful to make some additional remarks at this point to clarify the ideas 
involved. If a curve C C Ti is such that any two points on it (not too far apart) are in phase, 
then it is definitely horizontal: 




arg(^(s), i>(s')) = a(s') — a(s), any s' and s 
— — arg(^(s),^(s / )) =0 



arg(^(s), ip(s')) = separable in s' and s 




(3.4) 



C = {ip(s)} : (ip(s),ip(s')) = real positive 




C horizontal . 



(3.5) 
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The image C = n[C] is obviously a null phase curve in 1Z as eqn.(3.4) is obeyed with 
ct(s) = 0; therefore C being a lift of C is also a null phase curve in H. On the other hand, 
for a horizontal curve C H only "nearby points" are in phase: 



C (h) = |^,W( S )J = horizontal 



V fc) W,^ (fc) («))=o 



W( s ),^)( s + 5s)) ~ i + o(<5s) 2 
=^ arg(^W(s),^W(s + (5s)) =0(5s) 2 . (3.6) 

However two general points on may well be not in phase, as axg(ip(s),ip(s')) could be 
nonzero. Hence and its image 7r[C^] may not be null phase curves. For tt[CW] to be 
a null phase curve, in addition to being horizontal (a local property) must possess the 
global property that for general s and s' the inner product (ip( h \s), i(;^(s'fj is real positive. 
This is what is captured in the conditions (3.4). 

We can now generalise the result (2.9) and strengthen it as follows. Given n unit vectors 
■01, fa, ■ ■ ■ , ipn ^ H- with images pi,p 2 ,...,p n e 1Z : draw any null phase curves joining 
consecutive pairs of points p± to p 2 , pi to p 3 , . . . , p n to pi. (This can certainly be done since 
in any event free geodesies are available). Then by exactly the same arguments that lead to 
the connection (2.9) we obtain: 



-argA n (^ 1 ,0 2 ,...,0 n ). (3.7) 



n-sided figure in TZ with vertices pi, p 2 , . . . , p n 
and bounded by null phase curves 

It must be clear that this is the widest generalisation of the connection (2.9) that one can 
obtain. We see that we can replace each free geodesic belonging to a polygon in 1Z by any 
null phase curve, and the geometric phase remains the same, since the right hand side of 
eqn.(3.7) depends on the vertices alone. 
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IV. CONSTRAINED GEODESICS AS NULL PHASE CURVES 



We have seen that every free geodesic is a null phase curve, but the converse is generally 
not true. Nevertheless the former fact motivates the following: can we alter the definition 
of a free geodesic, based on minimising the length functional L[C] of eqn.(2.5), in a natural 
way to obtain other kinds of geodesies, and will they turn out to be null phase curves as 
well? 

The generalisation we explore is the following: instead of dealing with curves (of unit 
vectors) in the complete Hilbert and ray spaces Ji and 71, we restrict ourselves to some 
(smooth) submanifold M C 1Z and consider only curves C lying in M and connecting pairs 
of points in M. For such curves we minimise L[C] with respect to variations of C which 
stay within M. The resulting curves will naturally be called "constrained geodesies" , and 
the question is: do constrained geodesies in some cases turn out to be null phase curves? 

We emphasize that our question is not whether every null phase curve is a constrained 
geodesic lying in a suitably chosen submanifold M C 1Z, but rather whether the latter curves 
sometimes have the former property. The physically important examples presented in the 
next Section show that our question is indeed interesting. In this Section we set up the 
general framework to handle constrained geodesies in ray space. 

Given H and 1Z with dim H = dim 1Z + 1 in the real sense, we consider a submanifold 
M C 1Z of n (real) dimensions consisting of a (sufficiently smooth) family of unit rays, with 
(local) real independent and essential coordinates £ = (£ M ), ji = 1, 2, . . . , n: 

M = { P (0 e K | £ e W}cK. (4.1) 

(We do not indicate explicitly the domain in 3? n over which £ may vary). The inverse image 
of M in H will bring in an extra phase angle a, and is denoted by M.: 

M = 7r _1 [M] 

= {V(£;a) G H I 7r(^;a)) = p(0, m a) = e"W; 0)} . (4.2) 
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(Of course each a) is a unit vector, and a and £ M taken together are local coordinates 
for M. ). So in the real sense dim M. — and to avoid trivialities we must have l+n/2 < 
complex dimension of 7i. 

Now we consider a parametrised curve C <Z M G 1Z, obtained by making the n real 
variables £ M into functions of a real parameter s: 



To lift C to some C C M C H, some (smooth) choice of phase angle a(s) as a function of s 
must be made, and then we have: 



Using the definition (2.5) the length L[C] can be seen to involve only the partial derivatives 
of ■?/>(£; a) with respect to the the dependence on a being trivial and not contributing at 
all. Therefore we define: 



C = {p(£(s)), Sl <s<s 2 }cM. 



(4.3) 



C = {M/( s )=^(£( S );a( S ))}cM, 



7T[C] = C. 



(4.4) 



d 



il>{£;<x), /i = 1,2, ...,n; 



(4.5) 



Normalisation of a) to unity for all £ and a implies 



Re (V>(f;a),?v(f;a:)) = 0. 



(4.6) 



Now L[C] can be expressed as follows: 



L[C] 



fd8y/\\*(s)\\*-\(*(s)Ms))\ 2 




(4.7) 
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The parameter dependences of £ and a are as in eqn.(4.4). From the essentiality of ^ as 
coordinates for M, and the positivity of the metric on H, one easily obtains the following 
results: the nxn matrix (Juj^(^; a),u^(£; afj^j is hermitian positive definite and independent 
of a; and only its real part which is symmetric positive definite, enters L[C}. 

To obtain the differential equations for constrained geodesies we minimise L[C] with 
respect to variations in C that stay within M. This amounts to minimising L[C] in the 
final form given in eqn.(4.7), by making independent variations in the n real functions £ M (s); 
and the result is well-known from Riemannian geometry. After making a suitable choice of 
the parameter s (afline parametrisation), the differential equations for constrained geodesies 
become: 



(<T(0) = (<M0) _1 , 



9 P ,M = ™- (4.8) 



Here the T's are the familiar symmetric Christoffel symbols determined by the 'metric' tensor 
g^uiO- The remaining freedom in the choice of parameter s is only change in scale and shift 
of origin. It is a consequence of the differential equations above that 

9,MsW(s)C(s) = constant . (4.9) 

A general solution to eqn.(4.8) is uniquely determined by choices of initial values 
£ M (0), £ M (0). The resulting £ M (s) determine some constrained geodesic C C M C 71; and for 
any (smooth) choice of a(s) we get a lift C C M. C 7i which by definition is a constrained 
geodesic in H. The meaning of the 'conservation law' (4.9) in terms of Hilbert space vec- 
tors is interesting. In terms of the derivative of \l/(s) with respect to s, and its component 
orthogonal to \P(s), 

*(a) = ^lK£(a);a(a)) 
12 



*-L( S ) = *( S ) - *( S )(*(s),*(s)) 



(4.10) 



we have, 



9^(C(s))e(s)C(s) = constant 




^ f_L (s) ||= constant . 



(4.11) 



We can then, if we wish, adjust the scale of s so that \E'- L becomes a unit vector for all s. 



Having set up the basic formalism to determine constrained geodesies, in the next Section 
we look at some physically motivated examples to see whether they are sometimes null phase 
curves as well. 



We look at four examples to illustrate the use of constrained geodesies in the geometric 
phase context, and to show the distinction in general between them and null phase curves, 
(a) Subspaces of H 

Let Ho be a linear subspace of H (as a complex vector space), and denote by M. C Ho the 
subset of unit vectors in Ho- By projection we obtain the submanifold M = ir[M] C 1Z, 
with the real dimension of M equal to 2 x {(complex dimension of Ho) — 1}- In this case 
constrained geodesies in M happen to be free geodesies. For, given any two (non orthogonal) 
in phase unit vectors in A4, say ipi and ip 2 , the free geodesic connecting them, namely from 
eqn.(2.6) the curve C consisting of the vectors 



passes entirely through points of M. Hence its image ir[C] = C lies entirely within M; and 
being the free geodesic connecting n(if)i) to 7r(^ 2 ) it must be the constrained geodesic as 
well. In this case therefore we do not get anything new. 



V. APPLICATIONS 




sm s , 



(5.1) 
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Conversely we see that to have a situation where constrained geodesies are different from 
free ones, the submanifold M C 1Z must not arise from a subspace of H in the above manner. 
We now look at two such cases, of obvious physical importance, in which true generalisations 
of the original Bargmann invariant-geometric phase connection appear, 
(b) Single mode coherent states 

We consider the family of coherent states for a single degree of freedom, described by 
hermitian operators q, p or the nonhermitian combinations a, a) : 

1 1 

a=-j=(q + ip), a f = —(q-ip) , 

[q,p]=i, [a,tf] = l. (5.2) 

A general normalised coherent state is labelled by a complex number z and is generated by 
applying a unitary phase space displacement operator to the (Fock) vacuum state |0): 

\z) = exp(za^ — z*a) |0) 

= exp (^~2 Z * Z ~^ > 
a\z) = z\z) . (5.3) 

To conform to the notations of the previous Section, we introduce real parameters £1,^2, 
include a phase angle a, and express the above states in terms of q and p as follows (for ease 
in writing we use £1,2 rather than £ 1,2 ): 

iP(£;a) = e ia \z) 

= exp (ia + i (£ 2 <? - |0) 

= exp (ia- ^162) ex P(« 6?) ex P (~i 6p) |0) 

= exp (ia + ^£2) exp (-i ^p) exp(i £ 2 g) |0) . (5.4) 

(Note that, as in eqn.(4.2), ■?/>(£; a) is a vector in TL parametrised by £ and a, not a wave- 
function). These various equivalent forms facilitate further calculations. 
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The expectation values of q and p in these states are 

(^;a),p^;a))=&. (5.5) 

Now we compute the vectors a) and their projections a) orthogonal to a), 
as defined in eqn.(4.5): 



m 2 = d 2 il> = i(q- ip; 



2 

«i" = -i(p-6)^, «2" = «(9-6)^- ( 5 - 6 ) 

Here we used eqn.(5.5), and for simplicity omitted the arguments £, a in tpjU^uj^. The 
inner products among the uj; involve the fluctuations in q and p and the cross term. After 
easy calculations we find: 

(ui, ui) = (V, (p - 6) 2 V>) = (Ap) 2 = 1/2 , 
(< 4) = " (P " 6) (5 " 6) V>) = i/2 , 

(4, 4) = (V, (? - 6) 2 V>) = (Ag) 2 = 1/2 . (5.7) 
Therefore the induced metric tensor in the £1 — ^2 plane, defined in eqn.(4.7), is 

9 AO = \Kv , (5-8) 

namely it is the ordinary Euclidean metric on K 2 . Constrained geodesies in this case are 
just determined by straight lines in the £-plane, as all T's vanish: 

z(s) = z + zis, z ,i = ;^=(9o,i +iPo,i) : 

Ci(s)=5o + 5is, &(s)=p +p 1 s. (5.9) 

At the Hilbert space level, a constrained geodesic C constr .g C o. can be taken to be a curve 
within the family of coherent states 

15 



Cconstr.gco. = = \z + ZlS)} . (5.10) 

(Here we have omitted an s-dependent phase a(s)). Each vector ^f(s) along this curve is a 
(pure) coherent state, and can not be written as a linear combination of two fixed states as 
in eq.(2.6); so it is immediately clear that this is not a free geodesic at all. 

Now we examine whether this constrained geodesic is a null phase curve. We find, using 
the criterion (3.4): 

arg(^(s), *(s')) = arg(^ + Z!s\z + zis') 

= arg (exp {(z* Q + z{s) (z + z x sl)}) 
= arg (exp (zqZis' + z zls)) 

= (s' - s)lmz* oZl . (5.11) 

This is a separable function of s' and s, so we do have a null phase curve. We can go from 
the above C constr .geo. to a horizontal curve by adding a phase: 

CioL.geo. = {*» = exp {-i slm 4zi)*(s)} , (5.12) 

and then we find that any two points on this curve are in phase, as expected. 

The generalised connection (3.7) in this example now states: if \zi), \z2), ■ ■ ■ , \z n ) are any 
n pure coherent states given by choosing n points in the complex plane, and we join these 
points successively by straight lines in the complex plane, so that all along in Hilbert space 
we deal with individual coherent states and never with superpositions of them, we have: 



-argA n (ki), \z 2 ), \z n )) . (5.13) 



n-sided plane polygon with vertices 
at the coherent states z 1: z 2 , ■ ■ ■ , z n 

The case n = 3 leads to the area formula for the geometric phase for a triangle in the plane, 
a very familiar result ||15|| . From our point of view, the present example is a significant 
generalisation of the original connection (2.9). 

Going further, it is easy to convince oneself that in this example the most general null 
phase curve arises in the above manner; in other words, a given one-parameter family of 
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coherent states {|^(s))} obeys the separability condition (3.4) if and only if lmz(s) is a 
linear inhomogeneous expression in Rez(s), so z(s) describes a straight line in the complex 
plane as s varies. 

(c) Centred Gaussian pure states 

This example again deals with one canonical pair q, p. It is now more convenient to work 
with wavefunctions in the Schrodinger representation, and not with abstract ket vectors. 
The submanifold M. C 7i consists of normalised Gaussian wavefunctions parametrised by 
two real variables ^1,^2 and a phase angle a defined as follows: 

a; q) = (^j ex P + 1(6 + * 6)? 2 } , 

6 E (-00, 00), ^ e (0,oo), aG[0,27r). (5.14) 



Normalisability requires that £ 2 be strictly positive, so the combination ^ + i £ 2 is a variable 
point in the upper half complex plane. The wave functions ^(£5 a; g) are: 

d % 

; «; 9) = -^Kf ; «; 9) = « ^ 2 a; ?)> 

a; = a '^) = \ (^-v 2 + «; 9) • (5.15) 

It is clear that to obtain the components v,j; of w M orthogonal to if), and later to compute the 
inner products (u^,u^j, we need the expectation values of q 2 and q 4 in the state ip. These 
are (omitting for simplicity the arguments of ip): 



, p \ V2 °° 

(r.r/V) f^J j dqq 2 e-^ r — 



-00 



DC 



26' 



ty,9V) = f |J fdqq 4 e-^ 2 = -L . (5.16) 



Now the necessary inner products and projections are easily found: 

(^>«i) = 77-> (^,«2) = 0; 
4^2 
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2 
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U^) = (U 2 , M 2 ) = ^ f V, - V»j = • ( 5 - 17 ) 

From these results we obtain the induced metric over M = ir[Ai] C R, described in the 
upper half complex plane by the metric tensor 

9^(0 = Re a), u^(£; a)) = — . (5.18) 

This is the well-known form of the Lobachevskian metric in this model of Lobachevsky space 
Dropping the numerical factor 1/8 for simplicity, the line element in the upper half 



plane is given by 

ds 2 = ^ Ug + dZl) , (5.19) 

and we must find the corresponding geodesies. 

First we compute the nonvanishing T's. The inverse of (^(0) nas components 

9 1 \0 = 9 2 \0=e 2 , 9 U (O=0. (5.20) 
We easily find that the nonvanishing T's are 

r 1 1 2 (0 = r 2 22 (0 = -r 2 u (0 = - )- . (5.21) 

Using these in eqn.(4.8) we find the following ordinary differential equations to determine 
geodesies: 

6-1-66 = 0, (5.22a) 

6 + ^ Ui-H) =0. (5.22b) 

<52 v ' 

We can exploit the fact that these differential equations lead to the consequence 
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4fe 2 + 4 2 ) = constant, (5.23) 

the value of the constant depending on the particular geodesic. After elementary analysis 
we find that there are two families of geodesies: 

Type I: £1 = constant , £ 2 = ae bs , 

a > 0, s E 3?; (5.24a) 
Type II: ^ = c + Rcosf(s) , £ 2 = Rsmf(s) , 
f(s) =2tan~ 1 (ae 6s ), 

ceK, #>(), a>0, b>0, seS. (5.24b) 

These are both in affinely parametrised form. In Type II it is simpler to pass to a nonaffine 
angle type parameter s G (0,7r), and replace eqn. (5.24b) by: 

Type II: £i = c + R cos s, £ 2 = -Rsins, 

cGS, R>0, 0<s<n. (5.25) 

Type I geodesies are straight semi infinite lines parallel to the £ 2 axis. Type II geodesies are 
semicircles centered on the £i axis and lying above this axis. 

In each case we can now ask whether a constrained geodesic in M is a null phase curve. As 
in the previous example of coherent states, here too we emphasize that we are concerned with 
curves within the manifold of centred normalised Gaussian wavef unctions, and at no stage 
with linear combinations of such wavef unctions. We look at the two types of constrained 
geodesies in turn and find these results (after simple reparametrisations): 

Type I: = ^(6 = a > 6 = bs; a(s)) : 

arg(tf(s),tf(s')) = 0; (5.26a) 
Type II: ^(s) = ^(£1 = c + R cos s, £ 2 = R sin s; a(s)) : 

axg(*(s),*(s')) = J(s-s')- (5.26b) 
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(In both cases the choice of phase angle a(s) is irrelevant). So in both cases the criterion 
(3.4) is obeyed; and both types of curves in M arising from the two types of geodesies in 
the upper half £ plane are simultaneously constrained geodesies and null phase curves. 

The statement of the generalised connection (3.7) is clear, and for illustration we consider 
the case of just three vertices. Let A, A' and A" be any three points in the upper half complex 
plane; and for any choices of phases a consider the three normalised centered Gaussian states 
ip(A; a),ip(A'; a') and ip(A"; a"). Join A to A', A' to A" and A" to A by a geodesic of Type I 
or Type II as appropriate in each case. This can always be done, and we obtain a hyperbolic 
triangle. In M we obtain a 'triangle' with vertices p(A) = ir(tf}(A] a)) etc., and whose sides 
are constrained geodesies; and we can state: 



'triangle' in M with vertices 
p(A),p(A>),p(A") 
and sides as constrained geodesies 



arg A 3 (V(A; a),ip(A'; a'),^{A"; a")) . (5.27) 



An application of this result has been used elsewhere |T7| to show that the classical Gouy 
phase |18| in wave optics is related to a Bargmann invariant and hence is a geometric phase, 
(d) A subset of two-mode coherent states 

In the previous two examples we found that while constrained geodesies differed from 
free geodesies, they were nevertheless null phase curves and so led to important instances 
of eqn.(3.7). This is however fortuitous; the really important objects for our purposes are 
the null phase curves, and in a given situation constrained geodesies may well not be such 
curves. In our fourth and final example, dealing with a subset of states for a two-mode 
system, we will find that this is just what happens. However we will be able to completely 
determine all null phase curves directly, so that the generalisation (3.7) can be meaningfully 
stated. 

For a two- mode system with creation and annihilation operators df, dj obeying the stan- 



dard commutation relations 



dj , d k 



5jk , [dj,dk] 



dj, a k 



0, j,k = l,2 



(5.28) 
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the general coherent state is labelled by two independent complex numbers arranged as a 
column vector z = {z\, z 2 ) T : 

\z) = exp {^-\^ z + 2161 + ^4) |0>, 
a J -|z> = ^|«), j — 1,2. (5.29) 

Within this family of all normalised coherent states we now define a submanifold (of real 
dimension three including an overall phase), an "i? 2 - worth of states", by taking 6,4> to be 
spherical polar angles on a sphere S 2 and setting z\ and z 2 equal to the following: 

z l = cos6, z 2 = e^sin0, < 6 < n , < < 2tt . (5.30) 

Therefore the submanifold M. C is parametrised by 6>, and a phase a and we write: 

M = {ip(6,<f>;a) = e ia \cos9,e i(t> sin9) | < 6 < ix , < 0, a < 2tt} C H , 

(5.31) 

where the ket on the right is a particular two-mode coherent state with z^z = 1 

tp{9, 0; a) = exp (ia + a\ cos + c&e** sin - 1/2) |0) . (5.32) 

Omitting the arguments 6 1 , 0, a for simplicity, we easily find: 

d 

u e = —ip = (— sin 6 a\ + cos Q a£) ip , 
d 

u,p = —ip = ie 1 ^ sin 8 alip ; (5.33a) 

(ip,ue) = 0, (ip , u<p) = i sin 2 9 ; (5.33b) 

Ug = ug , = i sin 6 (e^al — sin 6j ip . (5.33c) 

Repeatedly exploiting the eigenvector relation (5.29) and its adjoint, we compute the inner 
products among the vectors in eq. (5.33c): 

(y,g,ufj = l, (ug ,u£) = i cos#sin#, (u^, u£) = sin 2 . (5.34) 
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Taking the real parts here, we see that the metric induced on M = ir[M] ~ S 2 in 1Z, 
parametrised by angles 9 and 0, is just the usual rotationally invariant one: 

goe(9,(f>) = 1 , 000 = 0, g^(9, 0) = sin 2 9 . (5.35) 

The corresponding constrained geodesies are therefore simply great circle arcs. The question 
is whether they lead to null phase curves in M and M. 

A general parametrised great circle arc on S 2 is traced out by an s-dependent unit vector 
n(s) with polar angles 9(s),<f>(s): 

n(s) = acoss + 6 sins = (sin#(s) cos0(s), sin#(s) sin0(s), cos#(s)) , 

a,beS 2 , a- 6 = 0. (5.36) 

The corresponding constrained geodesic C constr .g C o. C M. (omitting the phase a) is the curve 
of coherent states 

*(s) = \z 1 {s),z 2 {s)) , 
zi(s) = cos 9(s) = a 3 cos s + 6 3 sin s , 
z 2 (s) = e ilp{s) sin 9(s) = ( ai + i a 2 ) cos s + (&i + i 6 2 ) sin s . (5.37) 

To see whether this is a null phase curve we compute the phase of (^(s), fy(s')): 

arg (*(s), *(s')) = axg(zi(s), z 2 (s)Ms'), z 2 (s')> 

= arg{exp (zi(s)zi(s') + ^ 2 (s)*z 2 (s'))} 

= arg {exp ([(ai — « a 2 ) cos s + (6i — i 6 2 ) sin s] 

[(ai + i a 2 ) cos s' + (&i + i 6 2 ) sin s'])} 
= (a A 6) 3 sin(s' - s) . (5.38) 

Unless it vanishes, this is not a separable function of s' and s. We conclude that the geodesic 
(5.36) on S* 2 leads to a constrained geodesic C constr .g C o. C M. which is in general not a null 
phase curve. The only exception is when (a A b) 3 = 0, that is, the geodesic (5.36) on S 2 lies 
on a meridian of longitude, with a A b being a vector in the 1-2 plane. 
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On the other hand, in this example it is quite easy to explicitly find all null phase curves 
on M(and M)\ Let T = {n(s)} C S 2 be given, and let us consider the induced curve C r in 
M: 



C r = {^ r (s) = \n 3 (s),n 1 (s)+in 2 (s)) = 

exp(-^ +n 3 (s) a\ + (rii(s) + % n 2 (s))a|)|0)} . 



(5.39) 



We find that 



axg(tt r (s),*r(s')) = (ra(s) A n(s')) 3 ■ 



(5.40) 



This will be a separable function of s' and s if and only if, for some constants (5 and 7, we 
have 



The geometrical interpretation of this is that the projection of V on the 1-2 plane must be 
a straight line. In that case Cr is indeed a null phase curve in M., as we have 



which is separable in s' and s. One can easily see that each such T is a latitude circle arc 
on S 2 corresponding to (i.e., perpendicular to) some axis lying in the 1-2 plane; and given 
any two points on S 2 , we can always connect them by such a T. In other words, such T are 
intersections of S 2 with planes perpendicular to the 1-2 plane. When such a latitude circle 
arc is also a great circle arc, we recover the result of the previous paragraph. 

The upshot of this example is that here we have a nontrivial illustration of the difference 
between constrained geodesies and null phase curves. However, since we have been able to 
find all of the latter, and any two points in M. can be connected by some null phase curve, 
we have succeeded in providing a nontrivial two-mode example of the generalised connection 
(3.7), but not using constrained geodesies. 



n 2 (s) = (3 ni(s) +7- 



(5.41) 



arg(^ r (s), * r (s')) = 7M«) ~ n M)) , 



(5.42) 
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VI. RAY SPACE AND DIFFERENTIAL GEOMETRIC FORMULATIONS 



Very soon after the discovery of the geometric phase, the differential geometric expres- 
sions of its structure and significance were brought out |l9|,p|-0i|lTl , by relating it to anholon- 



omy and curvature in a suitable Hermitian line bundle on quantum mechanical ray space. 
In this section we provide a brief discussion of the properties and uses of the new concept 
of null phase curves at ray space level and also in the differential geometric language. Only 
necessary background material will be recalled, and derivations omitted. Since they may 
be useful for practical calculations, where possible local coordinate expressions of important 
differential geometric objects will be given. 

From the preceding sections it is evident that for our purposes it is important to deal 
with open null phase curves in general, since it is through them that the connection (3.7) 
of the Bargmann invariants to geometric phases is made. Their definition (3.4) in terms of 
Hilbert space lifts is quite simple. Nevertheless it is of interest to develop a direct ray space 
formulation; this can be done essentially via the Bargmann invariants themselves. From 
their definition (2.8), it is clear that any A2 is real nonnegative, while A n 's for n > 3 are in 
general complex. On the other hand it is also known that any A n for n > 4 can be written 
as the ratio of a suitable product of A 3 's and a suitable product of A 2 's: 

n n 

A n (Vi, Va, • • • M = II A 3(^i, V^-i, 1>j)l IT MMi-i)- (6-1) 

3=3 j=4 

In this sense the three-vertex Bargmann invariant A3 is the basic or primitive one as far 
as phases are concerned. (The basic cyclic invariance of A„(^i, ip2, • • • , ipn) is not manifest 
in eq.(6.1), but it is not lost either). Guided by these facts we give now a direct ray space 
characterisation of null phase curves. 

If C — {p(s)} C 1Z is a null phase curve and C^ 1 ' = {tp^ (s)} is a horizontal Hilbert space 
lift obeying eq.(3.4), we see immediately that for any choices of parameter values s, s', s", 

A 3 (^ h \s),^ h \s'),^ h \s")) = Tr{p(s)p(s')p(s")} = real and >0; (6.2) 

and so also for any n parameter values Si, s 2 , . . . , s n , from eqn.(6.1), 
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A n (^ h \ Sl ),...,^ h \s n )) = Tr{p( Sl )...p(s n )} = real and >0. 



(6.3) 



As a consequence, by differentiation with respect to s 2 , . . . , s n we have : 

dp(s 2 ) dp(s n )\ 



Tr {p( Sl )- 



real . 



(6.4) 



k ds 2 ds n J 

Now it is known that the geometric phase for any connected portion of any C can be 
expressed directly in terms of p(s) as follows, whether or not C is a null-phase curve: 

r*2 dp(s) 



(p g [p{si) to p(s 2 ) along C] = arg 



Tr < p(si)P ( exp J ds 



= art 



rS2 rs' n rs' 2 ( 

1 + 22 ds 'n / ds 'n-l ■ ■ ^i Tr S P(Sl 

n=1 Jsi J si J si [ 



ds 
dp(s[) 



(6.5) 



ds' n dsi 

where P is the ordering symbol placing later parameter values to the left of earlier ones. If 
eq.(6.2) holds on C (and so as a consequence eqs. (6.3,4) as well), we see that at every stage 
only real quantities are involved, the geometric phase in eq.(6.5) vanishes, and C is a null 
phase curve. This leads to the ray space characterisation of null phase curves we are seeking 



C = {p(s)} C 1Z is a null phase curve 

Tr{p(s) P (s')p(s")} 



real nonnegative, any s, s', s" . (6.6) 



Turning now to the specific differential geometric aspects, it is well known that while the 
dynamical phase V'dynJC] is an additive quantity, f g [C] does not have this property. On the 
manifold of unit vectors in Hilbert space TC, there is a one form A such that 



V?dyn[C] = [A. 

J C 



(6.7) 



However, referring to the projection n : 7i — > 71, A is not the pull-back via n* of any 
one-form on the space of unit rays; and (p g [0\ is not the integral along C of any one- form on 
1Z. In fact this lack of additivity can be expressed via the Bargmann invariant A 3 . If C\ 2 
connects p\ to p 2 in 1Z and C 23 connects p 2 to p 3 , than C i2 U C 23 runs from pi to p 3 and 



(Pg[C 12 U C 23 ] = (Pg[C 12 ] + (Pg[C 23 \ ~ B 3 (lp U tp 2 , ifo), 

^3(^1,^2,^3) = arg A 3 (^1,^2,^3)- 



(6.8) 
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More generally, for an (generally) open curve consisting of (n — 1) pieces Cu, C23, • • • C n -i,n 
joining p 1 to p 2 , p 2 to p 3 , ■ ■ • , p n _i to p n , we generalize eq. (6.8) to the following: 

n-l 

(p g [C 12 U C 23 U • ■ ■ U C n _i >n ] = J2 VgPjj+l) ~ B n(i>i, $2, ■ ■ ■ > Vv), 

i=i 

B n {lpl, 1p2,'--, ll>n) = ar g A n (^1, ^2, • ■ ■ , ^n) 

n 

= ^ J B 3 (Vi,^-i,^)- (6.9) 
i=3 

If we connect p n back to pi via C n> i to get a closed curve of n pieces, then we have the 
specific result: 

<f g [C 12 U C 23 U • • • U C„-l, n , UC„,i] = (p 9 [Ci 2 ] + ^[^23] H h Vs[Cn,l] - B^lf)!, l[) 2 , ■ ■ ■ , i>n)- 

(6.10) 

Compared to eq. (6.9), we have one extra ip g term on the right but the Bargmann phase 
term B n is the same. We see that the lack of additivity shown in all eqs. (6.8,6.9,6. 10) is due 
to the Bargmann pieces. There is however an exception to this general nonadditivity, which 
occurs in (6.8) when p 3 = pi and Ci 2 U C 23 is a closed loop. Then we find : 

<9(Ci 2 U C 23 ) = 0, p 3 = pi : 

<Pg[Cl2 U C21] = <Pg[Cl 2 ] + tPg[C 2 l] , 

i.e., <p g [C 12 ] = <p 9 [C 12 U C 21 ) - cp g [C 23 } . (6.11) 

In the past this result has been used Q to relate f g [C] for an open C to f g [C U C] for 
a closed C U C" by choosing C to be a free geodesic, for then <p g [C] = 0. Now we can 
generalise this process: if C is an open curve from p\ to p 2 in TZ, and C is any null phase 
curve from p 2 back to pi, we have the result 

<p 5 [open curve C] = (p g [closed loop C U C] . (6-12) 

This is the most general way in which an open curve geometric phase can be reduced to a 
closed loop geometric phase. More generally, comparing eqs. (6.9,6.10) valid for generally 
open and for a closed curve, we see that if the last piece C n> \ is a null phase curve we convert 
an open curve geometric phase to a closed loop geometric phase: 
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V 9 \Cv2 u c 23 u • • • U C n _i, n ] = ^[^12 u c 23 u • • • U C n _i, n U C n ,i] (6.13) 

At this point it is natural to express a closed loop geometric phase as a suitable "area 
integral" of a two-form, both at Hilbert and ray space levels. Whereas A is not the pull back 
of any one-form on TZ, we do have dA = it* uj, where a; is a symplectic (closed, nondegenerate) 
two- form on TZ. Then, if C is a closed loop in Ti, dC = 0, so that C = 11(C) is a closed loop 
in TZ, we have 

Vg [C] = f dA = [u, (6.14) 
Js Js 

where S and S = n(S) are two-dimensional surfaces in Ti and TZ respectively, with bound- 
aries C and C : 

dS = C, dS = C. (6.15) 

With the help of local coordinates on Ti and TZ we get explicit expressions for A, dA and 
uj. Around any point p £ TZ-, an d for some chosen ip G 7r _1 (po), we define an (open) 
neighbourhood iV C TZ by 

N = {p G TZ I Tr(pop) > 0} . (6.16) 

We can introduce real independent coordinates over iV as follows. Let {ip , ei, e 2 , • • • , e r , . . .} 
be an orthonormal basis for Ti. Then points in iV can be "labelled" in a one-to-one manner 
by vectors X G Ti orthogonal to ip and with norm less than unity: 

IIX(A7)H 2 = ^E(^ 2 + 7, 2 )<1: 

7) = x(/9, 7) + y/l-lWnWik , 
p e N <^> p = '0(/5,7)W,7) t , forsome/5,7. (6.17) 

Thus the real independent /3's and 7's, subject to the inequality above, are local coordinates 
for N. They can be extended to get local coordinates for 7r _1 (A^) C Ti by including a phase 
angle a : 
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g(v) = 1 + o i — r^F~ + o Jr iv J, 



i> E tt~ 1 (N) i> = i/;(a; (3, 7) = e ia ip(/3, 7) , < a < 2n . (6.18) 

In these local coordinates over N and 7r -1 (./V) we have the expressions 

1 

A = da + - ^2(-f r dp r - prdjr) , 
z r 

dA = ^2 d>lr A d(3 r , 

r 

= d"fr A d/3 r . (6.19) 

r 

The closure and nondegeneracy of u are manifest, so it is a symplectic two-form on TZ; and 
the coordinates f3, 7 realize the local Darboux or canonical structure for it. On the other 
hand, in these "symplectic" coordinates the Fubini-Study metric is a bit involved. If we 
combine the /3's and 7's into a single column vector r\ = ( (5\ P2 ■ ■ ■ 7i 72 • • • ) T , then the 
length functional L[C] of eq.(2.5) assumes the following local form: 

L[C\ = J ds^jrf g{rf)r] , 

1 7]1] T 1 

2 1 - ±?? T 77 + 2' 
1\ ' 

J = I , r] T r] < 2. (6.20) 

-1 0/ 

The symplectic matrix J plays a role in this expression for the metric tensor matrix g(i]). 
This matrix g(i]) is checked to be real symmetric positive definite, as one eigenvalue is 
(1 — \rj T rf)' 1 (eigenvector rf), another eigenvalue is (1 — \rj T rj) (eigenvector Jrj), and the 
remaining eigenvalues are all unity. We appreciate that for considerations of geometric 
phases and null phase curves this kind of local description is really appropriate, while free 
geodesies appear unavoidably complicated. 

We also notice that, in case TC is finite dimensional and the real dimension of the space 1Z 
of unit rays is 2n, the symplectic two-form u of eq.(6.19) is invariant under the linear matrix 
group Sp(2n, R) acting on the local coordinates (3, 7. On the other hand, the integrand of 
the length functional L[C] in eq.(6.20) possesses invariance only under Sp(2n, R)DSO(2n) ~ 
U(n), which is just the group of changes in the choice of the vectors {e r } which together 
with ip make up an orthonormal basis for H. 
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Returning now to the discussions in Sections III and IV, we can bring in submanifolds 
M C 1Z, M. = 7r _1 (M) C H, with local coordinates a as indicated in eqs. (4.1,2). Let 
iM '■ M. ^> 7i and %m '■ M <^-> 1Z be the corresponding identification maps. Straightforward 
calculations show that the pull-backs of A, dA, uj in eq.(6.15) to M. and M are locally given 
(with mild abuse of notation) by: 

i* M A = da + Im a), u^; a)) d^ , 
i* M dA = i* M u = lm(u^;a),u^a)) Ad? 

= Im «),«„(£; a)) d^ A d? . (6.21) 

We see, as is well known, that while the real symmetric part of the hermitian matrix 
(i u ^i u v)) determines the metric, eq.(4.7), the imaginary antisymmetric part of the same 
matrix is relevant for symplectic structure and geometric phase, reinforcing the link between 
the latter two. (In case M — 1Z and M. = H, the £ M 's become the /3's and 7's of eqs. (6.17), 
and we immediately recover the expressions (6.19)). For our present purposes the follow- 
ing comments are pertinent: While uj is closed and nondegenerate, ImOU is closed but may 
well be degenerate. An extreme case is when M is an isotropic submanifold in 1Z, for then 
i* M uo = 0. Such a situation can easily arise if, for example, M is described by a family of 
real Schrodinger wavefunctions i/)(£;q)- (A Lagrangian submanifold in 1Z is a particular case 
of an isotropic submanifold when the dimension is maximal, namely half the real dimension 
of 1Z). One may expect that if M is isotropic and C C M, then C is a null phase curve. 
However this need not always be so, and the situation is as follows. For a general open 
curve C12 from pi to p 2 in a general submanifold M, if we can find a null phase curve C 2 i 
from p 2 to p\ also lying in M, then C\i U C21 is a closed loop; if ni(M) = 0, we can find 
a two-dimensional surface S G M having C12 U C21 as boundary. Then from eqn.(6.12) we 
obtain under these circumstances 

<Pg[Cl 2 ] =<P g [Ci 2 UC 2 i] = ! %uU. (6.22) 

JseM 

Here as stated above, we had to choose C 2 i to be a null phase curve. (In case p 2 = p\ and 

29 



C12 is already a closed loop, there is no need for any C 2 \\ it can be chosen to be trivial!) If 
however M is an isotropic submanifold, i.e. i^oj = (and assuming also ni(M) = 0), we 
can extract some very interesting consequences for geometric phases, though it falls short of 
the vanishing of <p g [C\ for every C C M. We have the chain of implications 

i^u = iljoo = 0, any two dimesional S C M 
J s 

<P g [Ci2 U C 2 i] = 0, anypi, p 2 , C 12 , C 21 in M 
^ ^fg[C 12 ] unchanged under any continuous deformation 

of C12 leaving the end points pi, p 2 fixed (6.23) 

Thus, within an isotropic submanifold, the geometric phase for a general curve depends on 
the two end points alone. In case the curve chosen is closed, it can be continuously shrunk 
to a point ( since, tti(M) = 0) and then its geometric phase vanishes. One can thus say in 
summary: 

McK,i£u; = 0, tti(M) = 0, CcM: 
ipg[C] = if dC = 0; 

if g [C] = function of DC alone, if dC ^ 0. (6.24) 

The main conclusion is that general open curves in an isotropic submanifold need not be 
null phase curves; but geometric phases are invariant under continuous changes of their 
arguments leaving the end points unchanged. Perhaps this is not too surprising after all, 
since the isotropic property is a two-form condition. 



VII. CONCLUDING REMARKS 

We have shown that the familiar connection between the Bargmann invariants and ge- 
ometric phases in quantum mechanics, based on the properties of free geodesies in ray and 
Hilbert spaces, can be generalised to a very significant extent. In fact we have shown that 
our generalisation is the widest possible one. The essential new concept is that of null 
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phase curves in Hilbert and ray spaces - the replacement of free geodesies by such curves 
leads to our generalisation. We have seen through examples that this wider connection be- 
tween Bargmann invariants and geometric phases is just what is needed in several physically 
relevant situations. 

Motivated by the fact that free geodesies are always null phase curves, we have defined 
the concept of constrained geodesies and posed the problem of determining when these may 
be null phase curves. We have presented two examples when this is indeed so, and one 
where they are not the same. This reemphasizes the fact that constrained geodesies and 
null phase curves are in principle different objects, and sharpens the question of finding 
useful characterizations of the former which may ensure the latter property for them. This 
is sure to shed more light on the general questions raised in this paper, and we plan to return 
to them elsewhere. 
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